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A general discussion of graphene in external electromagnetic field is provided. In general, the formulation 
is not Lorentz invariant because of Zeeman energy. But it can be restored approxiamtely in the case of strong 
' magnetic field, the condition when quantum Hall effect is observed. Besides obtaining the well-known Hall 

04 . conductance + 1/2), we also provide an explanation of the newly observed Hall conductance 

P5 I for L = 0, 1. These are part of the sequence of Hall conductance which depends on the filling of 

^ . Zeeman levels. The energy levels are obtained for general orthogonal constant electric and magnetic field. 

The second order Dirac equation is derived in applied monochromatic electromagnetic wave, and the major 
I difference between graphene system and convetional Dirac electrons is pointed out. 
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QJ , I. INTRODUCTION 

Graphene has been one of the major foci in physics because of its simple lattice structure and linear dispersion relation 
Ci • near the Fermi level [l]-[3] when only nearest-neighbor hopping is taken into account. It has become a new testbed not only 
for condensed matter physics, but also for quantum field theory and mathematical physics [4] -[5]. The physical properties 
of graphene in external field such as quantum Hall effect (both integer and fractional [6]-[ll]), spin quantum Hall effect [12], 
transport theory [13] [14], superconducting [15], and magnetic confinement[16] are under intensive study. Effects of next-nearest- 
neighbor hopping have also been studied [17]. It is widely recognized that the integer quantum Hall conductance is a^y = 

O , + i) where n = 0,l,---[l]. Although disorder and 4-fold symmetry breaking may be utilized to explain the newly found 

quantum Hall structures v — 0, ±1, ±4[18] -[22], a unified explanation is still called for and the simple structure of axy — 

. ^("- + 5) still deserves a simple and fundamental explanation. Not surprisingly, quantum Hall plateaus 1^ = ±2, ±6, ±10 
J> ' in graphenes can be explained by using the Landau levels of spinless particle in external magnetic field and the 1h-2 Lorentz 
invariance of the massless Dirac Hamiltonian [23]. Because simple utilisation of Landau levels does not consider the Zeeman 
energy, which is not negligible compared with low-lying Landau levels, Zeeman energy might explain some of the newly found 
plateaus. Yet, Zeeman energy is not Lorentz invariant. Therefore, how to use the well-known Landau levels and the Lorentz 
transformation property of Dirac Hamiltonian remains an issue. In this paper, we discuss graphene in general external fields. 
Because of the Zeeman energy, the system is not SO(l,2) invariant in general. But when the magnetic field is strong enough, the 
system does enjoy the invariance approximately, which can be utilized to relate the physics of one experimental configuration 
to that of another. We obtain the energy levels of the system for the case of constant magnetic fields, and for constant electric 
and magnetic field. Finally, we derive the second-order Dirac equation in monochromatic electromagnetic waves and point out 
the difference bewteen graphene system and conventional relativistic electrons. The rest of the paper is organized as follows: 
in Section II, we present the Hamiltonian in general applied field in the nearest-neighbor hopping approximation and discuss 
jj] condition of SO(l,2) invariance of the corresponding Lagrangian. In Section III, we discuss the energy levels of the Hamiltonian 
Ci in applied perpendicular magnetic field. We dicsuss in Section IV the energy levels in applied orthogonal magnetic and electric 
fields and derive the quantum Hall conductance in the case of strong magnetic field using the approximate SO(l,2) invariance. 
A discussion of Dirac equation in applied monochromatic electromagnetic wave is also provided using the method of VoUcow in 
this section. The last section is a brief summary. 
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II. HAMILTONIAN IN APPLIED FIELD 

The direct lattice of graphene is a superposition of two interpenetrated triangular lattices A^, A^. The generators of lattice 
Aa are [24] ai = \/3a(|, ^j^), and = V^a{^, where a « 1.42 A is the carbon-carbon distance. The vectors 

Si = a(0, —1), S2 = ^) and S3 = o,{^^, i) connect each site in the lattice A^i to its nearest neighbor sites in the lattice 

Ab. Unlike regular electron spin, the pseudospin in graphene represents the two sublattices and there is no magnetic moment 
associated hence does not couple directiy to magnetic field[6]. The tight binding Hamiltonian can then be written 

where a is pseudo-spin index and t is the uniform hopping constant. In the presence of applied magnetic field B and electric 
fields E = — V(p, Zeeman energy and Coulomb energy should be included. 

Hz = Mb[E B-A\r)TA{r)+ ^ B ■ B\r)rB{r)] (2) 

Hc = q[Yl ^{r)A^{r)A{r) + ^ ip{r)B'< {r)B{r)] (3) 

, fig = is the Bohr magneton. In momentum space, with ttfe = ^= J2 e~*'°'^A(r), 6fe = J2 e~*'°"^B(r), where 

A/a is the number of lattice points of sublattice A^ (or A^), Hq + Hz reads in fe-space 

H = -t Y^ifkalkbak + fkblka^k) + mB ■ ^{alrak + blrhk) (4) 

k k 

with fk = Ei=i e'**'**. t w 3.033 eV. 

/fe = e-^'^"" + 26^'="" cos(^fc,a) (5) 

The lattice Hamitonian Hq vanishes at the six comers of the first BriUouin zone. Among these, only two are inequivalent, and 
can be chosen as 

Let k = K± + p/h. In the vicinity of K+, we have 

3 3 3 
fk = --^Pxa/h - -ipya/h + -a'^ipx - Wyf/fi'^ H (7) 

In the vicinity of K_ , we have 

3 3 3 

fk = -^Pxa/h - -jiPya/h+ -a^{px + ipyf/h^ H (8) 

After defining two-component spinors as V'fe = (ofe, &fe)^, V'l = ^fc)> V'feai = o-ka, V'feo-2 =bka,<^ = ±1- we have 

^ = E ( f* '^n ) ^'"^ + • ^(V'LlO-«/3V'fc/31 + V'L2'^"/3V'fc^2) + He 



(9) 



fecr 

When Ho is linearized around these two points , one obtains 



Hoik=K^+p = v,(^^^^.p^ ^-y^^) =^,a.p (10) 



where = ^,cx = {ax, -o-y). and 



Px ~iPy\ _ 



-v^cr p 



<^ = (f X5 cTy) • In the case B = (0, 0, B), we have 



(11) 



(12) 



Or,letVfe+ = a= { 



, V'fe- = % = 



, we have 



(13) 



(14) 



Close to K+, we have 

p p 

There are two inequivalent representations of the 7-matrices in three dimensions:7'' = (fa, i<J2, icri) and 7'' = 
{—(^3, —i(T2, —icri)- We choose the first. 

H = Y1 fp(^F7 • P + Mb^7°)?p + ^p(^f7 • P - l^uBj°)Vp + He (15) 
p p 

where ^ = ^^7". Incorporated U{1) gauge invariance, the Hamiltonian reads 

H = J d^x^{x)[hv^j-{-iV-^A) + n^Bj° + q^°^]C{x)+ J d^xfj{x)[hv^-f ■ {-iV - ^A) - fi^B^° + q^°ip]r]{x) (16) 

Let a;'* = {x°,x) = {v^^t, x, y),A'' = A). Denoting D^^df, + ij-A^, = if/v^, we have 

C = ^hhdtC -n= C{x){ihp + gBj°)^{x) + fj{x){ihp - gBf)r]{x) (17) 

where g = 11^ /v^. The U(l) gauge invariance is preserved but the Lorentz SO(l,2) invariance is broken by the Zeeman term. 
Since 



For \E\ « V \B\ , we can write 



and in this case 



-e^^-d^A^^r = B7° + —Exj^ - —Ey^' 



^{x){ihp - ge^'''^d^A,^r)^{x) + 7){x){ihp + ge'^'^^d^A^^rHx) 



which shows 1+2 Lorentz invariance. The current is 

SL 



(18) 



(19) 



(20) 



(21) 
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ni. CONSTANT PERPENDICULAR MAGNETIC FIELD 



In this case cp = and we want to calculate the grand-canonical partition function. Z = Tre ,K = H — fiN. The Dirac 
equation is 



hv^a ■ (-iV - -A) + jJL^B — jjL 



which turns into second order 



{-DiDi + -Bas)^ € 

Using standard Landau levels (assuming qB = \qB\), we have 

{K-tx^B + iif = hvl2qB{£ + 1/2 - s,) 

So 



The partition function is 



Ke^s, =±hv^\j2^B{£+ 1/2 -s^) + H^B-ij. 



where = is the Landau degeneracy per unit area. For B = IT, we have « 2.4 x 



-r = A,^^ln[l + e- 



For = 0, we have the energy levels 



Corresponding to field r], we have 



Ke,s^ = ±hv^^j 2 j^B{e+ 1/2 - s,)+ hbB 



Ke,s^ = ±hv^^/ 2 j:B{i+ 1/2- s,) - (ibB 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



So the energy levels are symmetric under Therefore if the graphene is undoped, the Fermi level is still at /z = 0. 

Suppose fisB > 0. For ^-field, levels are {[A] represents integer part of A) 



fiBB - h\vj^2fm., £ = 0,1, 



[^], positive 



1/2 :< ^^B - h\v,\^^m., £=[^] + l,... ,oo, negative 



(30) 



IJ.^B + h\v^\^/2fB£, £ = 1,2, ••• ,oo,positive 



levels are 

s^ = -l/2: { 
For ry-field, -?^£,i/2 levels are 

= 1/2 



ti^B-h\v^y2^B{i+l), ^ = 0,1,--- ,[^]-l, positive 

2 ^ 

^i^B-h\v^\^2^B{e+l), ^= [^],--- ,oo,negative 
I H^B + h\Vp\^/2^B{i + 1), ^ = 0, 1,2,- --oo, positive 



-H^B + h\v^\./2fm, £ = 0,1,--- ,[||^],negativcc 



(31) 



-^l^B + h\v^\^^m, £= [2^] + l,--- ,c^, positive 
-ji^B — h\v^\^j2^Bl, £ = 1,2, - ■ ■ ,00, negative 



(32) 
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-fi^B + h\v,\^2'iB{e+l), £ = 0,1,--- ,[^]-l, negative 



-^I^B + h\v,\^2!^B{e+l), £=[^],... ,c^, positive (33) 
-^l^B - h\v^ I ^2|5(£+ 1), ^ = 0, 1, 2, • • • oo, negative 

So negative levels comes from both ^ and rj fields. 

^:Ke,i/2:f^^B-h\v,\^l2^i,i=[^^] + l,--- ,<x^ (34) 



^:Ke,-^/2:fi^B-h\v,\^l2^B{£+l),£=[^^],... ,cx> (35) 



r?:i^,i/2:<! "/^b^ + ?lKlv/2pZ, i = 0,1,. . . ,[^] 
1 -Mb-B - %plv/2p?, ^ = 1, 2, • • • , 00 



r/:if, 1/2: <! -feg + %pl V2f 5(^ + 1) , e = 0,l,...,[^]-l 
\ -ii^B-h\v^\^2^B{i + l), ^ = 0,l,2,---oo " 



(36) 



(37) 



For t = 3.033 eV, = ^ « 106(m/s). Defining Bq := then « 1.1 x lO^T. For lab field B ~ lOT, hence, [^] = 0. 

Thus for ^-field, levels are 

{^^B 1 = 0, positive 

p,^B-'h\v^\s j2^Bi , I,--- ,(X),negative (38) 
[i^B + h\v^\-^2^Bl, ^= 1,2,- •• ,00, positive 

levels are 



„ - _i /o. / /^B-S-^|^p| V2f-B(£ + l) , ^ = 0,-- - ,00, negative 

' ■ \ /i3S + ;i|t;p|^2f 5(^ + 1), ^ = 0,1,2, •••00, positive ^ ' 

For ry-field, K 1^1/2 levels are 

{—figB, ^ = 0, negativee 

-At^B + h\v^\ ^2^Bi , £ = 1, • • • , 00, positive (40) 
-/UbB - h\v^\-^2\Bi, e = l,2,--- ,oo,negative 

Ke,-i/2 levels are 



, _ ,^J-ii^B + h\v,\ ^2j-B{i + l} , £ = 0,--- ,00, positive 

' ■ \ -/i3B-?l|t;p|V2fB(^+l), ^ = 0,1,2,- --oo, negative ^ 

Note = 1/2 belongs to sublattice Aa and Sz = — 1/2 belongs to sublattice A b ■ Consider negative levels 
For ^-field, levels are 

Sz = 1/2 : { n^B - h\v^\,/2lm, £ = 1, • • • , 00, negative (42) 

Ke^-i/2 levels are 

Sz = -1/2 : { fi^B - h\v^y2^B{e+l), ^ = 0, • • • , 00, negative (43) 
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For Tj-field, levels are 



Sz = 1/2 : 



— /XgiJ, £ = 0, negatives 

-/Xjji? — fi|Up|^2^i?^, € = 1,2, ••• ,00, negative 



levels are 



So levels 



and levels 



s^ = -l/2: { -f,^B-h\v^y2^B{i + l), € = 0, 1, 2, ••• oo, negative 



(44) 



(45) 



/x^S-;i|t;p|^/2^Bf,€= I,--- ,oo 



-fj,^B - h\Vp\^2-^Bi,i= 1,2,--- ,00 

are doubly degenerate(apart from the Landau degeneracy). We can not have an infinite number of negative levels fiUed-up 
since a cut-off is necessary. Suppose among the levels ji^B — Sjwp | y^2|B?, I = L"^ ■ ■ ■ ,L\ are filled, and among the levels 

-fx^B - fi\v^ I ^/^Be, f = L^, • • • are filled. If -n^B is also filled, then at T = OK 



N = AJ2(4 - L« + 1) + 2(4 - + 1) + 1) = 2AJLi + ^2 + ^) 

where Li = L'j^ - + 1, L2 = ^ if + 1- If ^Mb-S is not filled, i.e., /i < -n^B, then 

TV = A, (2(4 - 4 + 1) + 2(4 - 4 + 1)) = 2AJL1 + L2) 
The physics near K- makes equal contribution. The magnetization is 

M = /z3(iV^-iV,) 

where 



(46) 

(47) 
(48) 
(49) 



The effective potential for ^-fields should be regularized. The sum in the second term in the fowUowing formula actually extends 
to only Li instead of to infinity. Hence 



= In 



00 1^1 



e=i 



(50) 



With {N) = ^-^^^ = -^"^i^, we have 

1 



Similarly 



(A^)r, = A, 



2E 



1 



+2E 



1 



At zero temperature, /x = 0, and we have 

(iV) = (Ar)j + (iV)^ = AJl + 2(Li + L2)] 

and the manetization 

M = A^/i3(2Li-2L2-l) 

Since Li = L2, L2 + 1, we have M = ±A^ij,^ . 



(51) 

(52) 

(53) 
(54) 
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IV. PERPENDICULAR MAGNETIC AND ELECTRIC FIELDS 
A. Constant fields: \E\ « Vp\B\ 

In this case, we can use the 2+1 Lorentz invariance and the results of previous section. Suppose system S' is moving at 
velocity v relative to system S. For v = {v, 0), the SO(l,2) transformation is 

A'^ = ''A, (55) 
/ 7 -^7 0\ 

A^'^= -/37 7 (56) 

V i; 

where^=^,7 = ^. 

/ E,/v, Ey/vA 

Fi,^, = -EJv^ -B, (57) 

\-Ey/v, j 

The electromagnetic fields are related by 

= E\\ , E'j_=-f{E + vx B)± (58) 

Bfl = B|| , = 7(5 - ;J X E)i_ (59) 

Suppose we have inY, : E = 0,B = {0,0,B), then in S'. 

=0, E'j_=jvxB (60) 

B[l=B||, B'^=jB^ (61) 

We want = (0, E'y,0),B' = (0, 0, -B'). Hence we can choose v = {v, 0). In this case. 

E' = -jvB = 75 (62) 



So we have 



Soin E 



Note that 



The Hall conductance is 



if — /Ugi? is filled or 



« = (63) 



f'=if,r = -70f=i^f (64) 



f = q{v,p,j) (65) 



+ (66) 



<y = ^/ = 'j^L (67) 
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if — /Ug-B is not filled, where L = Li + L2. The physics near K- makes equal contribution, hence 



V B' 



4g 



-/i^B is filled 



(68) 



[ \L, is not filled 

The aboveSO(l,2)transformation breaks down when t; > v^, which means when \Ey\ > v^\B'^\. The magnetic moment current 

js = Mb (P€ - Pri)^k - jri) (69) 

we then have 



qE' 

if = (-7/3KM = M3^ 



(70) 



So far the sequence of Hall conductance axy = + 5) has been observed, we here predict the existence of the sequence 

<^xy = aiid the recently observed L = 0, 4 are just part of this sequence[22]. According to above analysis, the filling of 
Zeeman levels makes the difference. 

B. General orthogonal constant B and E 

When the condition \E\ « v^\B\ does not hold, we need to solve the Dirac equation directly. Let A = {—yB, 0), = yE. 
The second order Dirac equation reads 



DiDi + -Ba +-^<^' ^2 



Denoting e=|^,K=^, 



(71) 



AA + + iea^ - {ey - 



We can diagonalize the part + iea^ by Jordan decomposition (a similarity transformation). 



where 



(72) 



(73) 



(74) 



Let 



J 



^47r^A^ - £^ 



Ji 
J2 



(75) 



(76) 



Then 



AA + J-{ey- nf 



* = 



(77) 



r,2 «2 2igB q^B^ 2,7/ x2 

5y - 5^ ^^^^ + + J - (ey - k) 



* = 



(78) 



- 5^ - AnA^yid^ + 47r^A^j/" + J - (ey - k 



^' = 



(79) 
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Let 



Denoting = Px/h. 



Let 



Let 



then 



dl + kl+ AnA^ykx + An^Ay + J - {sy - 

b = AnA^kx + 2eK 
{-d^ + ay^ + by + c)h=0 



{-dl + au^ + c-—)h=0 
4a 



/ = 



c = kl-K'^ + Ji 



(80) 



(81) 



(82) 



(83) 



(84) 



(85) 



When a > 0, this is a hannonic problem while for a < 0, Ji imaginary. Similar issue exists in [25], thereof /3 can be larger 
than 1, hence 7 is imaginary. The critical value of B is Be = E/v^. For B < the solution should be parabolic cylinder 
functions, as in the case of [26]. 

For positive a, the problem is harmonic oscillator with m = 7i^/2, w = h. Hence we have 



We have 



For /2, = Va, we have 



for /i, Ji = —y/a, we have 



K = 



b'^ - 4ac = Aahuj{n + 1/2) 



27rA,. 



K = 



-gfc^ifc v/2a3/2(n+l) 
27rA,. 



^= 2^^ 

Now the energy levels depend on kx and Landau degeneracy is removed partially. 

~ekx ± V2oNhi 



E±{kx,n) = \x^B + hv^ 



277 A, 



(86) 



(87) 



(88) 



(89) 



(90) 



C. External monochromatic electromagnetic wave 



Consider an external vector potential 



which represents a monochromatic plane wave and satisfies Lorentz condition 



(91) 



a^A^ = 



(92) 
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. k = ^. Notice that the plane wave is of speed c, not Up . Hence 



kx Eq = cBt 







(93) 



theniSo >> WpBo- fc^ = i^,k),k = \k\k,x'^ = {vpt,x). Note fc'^/e^ = ^ ~ I'^P 0. This is a major difference between 

F F 

graphene system and real Dirac sysytem in appUed electromagnetic wave. We ignore the Zeeman term at first. 



ihp^ = 

As in [28], denoting (j) = k ■ r — uit and = (</>). Hence 



Or 



The filed tensor 



The second-order equation is 



Using 



and 



we have 



k-A:= k^.A^' = 



[7^71M'. = 2(7•fc)(7•^') 



2ilA>^d^ 

h 



^A^ + il(^.k){^-A) 



We seek a solution of this equation, a la the original Volkow ansatz [27] [28] 
then {A' = iA) 

2q 



k^F" - p^F - 2ik ■ pF' + -^A ■ pF 

n 



^A'F-^{j-k){j-A)F = 



(94) 



(95) 



(96) 



(97) 



(98) 



(99) 



(100) 



(101) 



(102) 



(103) 



We can impose (since we can always add to p a multiple of k to meet this condition yet the the functional form of ^ remains.) 

= or some other condition to simplify the equation (103), but it will always be a second-order equation. The reason is that 
here k^k^ ^ 0, as distinguishes graphene system from the conventional relativistic electrons. 



V. SUMMARY 



To summarize, we discussed the Hamiltonian and energy levels of graphene in general constant external electric and magnetic 
fields. The systems is not SO(l,2) Lorentz invariant when Zeeman energy is taken into account. But when the magnetic field 

is strong enough, SO(l,2) Lorentz invariance is well preseved. Employing the symmetry, we predicted a sequence a^y = 

and explain the recently observed Hall conductance axy = L = 0, 4, which is an indication that the Zeeman levels are not 
filled at zero temperature. The second-order Dirac equation is derived when the applied field is a monochromatic electromagnetic 
wave and the difference between graphene sysytem and standard relativistic electrons in this case is revealed. 
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